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ABSTRACT 

We summarize the analysis of a uniform set of both previously-known and newly-discovered scattered-light 
echoes, detected within 30" of SN 1987 A in ten years of optical imaging, and with which we have constructed 
the most complete three-dimensional model of the progenitor's circumstellar environment. Surrounding the 
SN is a richly-structured bipolar nebula. An outer, double-lobed "peanut," which we believe is the contact 
discontinuity between the red supergiant and main sequence winds, is a prolate shell extending 28 ly along the 
poles and 1 1 ly near the equator. Napoleon's Hat, previously believed to be an independent structure, is the 
waist of this peanut, which is pinched to a radius of 6 ly. Interior, the innermost circumstellar material lies along 
a cylindrical hourglass, 1 ly in radius and 4 ly long, which connects to the peanut by a thick equatorial disk. The 
nebulae are inclined 41° south and 8° east of the line of sight, slightly elliptical in cross section, and marginally 
offset west of the SN. The 3-D geometry of the three circumstellar rings is studied, suggesting the northern 
and southern rings are located 1.3 and 1.0 ly from the SN, while the equatorial ring is elliptical (b/a < 0.98), 
and spatially offset in the same direction as the hourglass. Dust-scattering models of the observed echo fluxes 
suggest that between the hourglass and bipolar lobes: the gas density drops from 1-3 cm" 3 to > 0.03 cm -3 ; the 
maximum dust-grain size increases from ~ 0.2/im to 2/j,m; and the silicatexarbonaceous dust ratio decreases. 
The nebulae have a total mass of ~ 1.7M©, yielding a red-supergiant mass loss around 5 x 10~ 6 M Q yr" 1 . 
We compare these results to current formation models, and find that no model has successfully reproduced this 
system. However, our results suggest a heuristic evolutionary sequence in which the progenitor evolves through 
two "blue-loops," perhaps accompanied by a close binary companion. 

Subject headings: circumstellar matter — dust — scattering — stars: mass loss — supernovae:individual (SN 
1987A) — techniques: image processing 



1. INTRODUCTION 

Supernova 1987 A, the first naked-eye supernova (SN) in 
four centuries, was first discovered on 1987 Feb. 24.23 at 
the Las Campanas Observatory (Shelton 1987), and was 
later identified with the Large M agellanic Cloud (LMC) pro- 
genitor, Sk -69° 202 <Sandulealdli%9T) . a B3 I supergiant 
jRousseau et alJfT978l) . From its spectral type and distance, 
we infer a luminosity of (3 - 6) x 10 38 erg s" 1 , a surface tem- 
perature 15000-18000 K, and a radius (2-4) x 10 12 cm, 
which imp ly that the star exploded as a blue supergiant (or 
BSG; IWooslev etaDll987l) . Spectra taken shortly after the 
explosion was discovered revealed Balmer lines, indicating 
a Type II SN, but expanding at nearly 0.1c. A number of 
anomalies, such as this high expansion velocity, a weak ini- 
tial outburst, a very hot UV spectrum at early times, and sub- 
sequent rapid color evolution, indicated a compact progenitor 
envelope, typical of a BSG, rather than the expected red super- 
giant (RSG) progenitor. Heavy mass loss could have driven 
Sk -69° 202 to the BSG stage (mass loss being reasonable 
given the large number of Wolf-Rayet stars in the LMC), al- 
though evidence for a substantial envelope (the long slow rise 
of the optical light source, the late appearance of near-IR H 

1 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, 
MD 21218; sugerman@stsci.edu 

2 Department of Astronomy, Columbia University, New York, NY 10027; 
arlin@astro. columbia.edu 

3 Guest Observer, Cerro-Tololo Inter-American Observatory 

4 Las Campanas Observatory, Carnegie Observatories, Casilla 601, La Ser- 
ena, Chile; kunkel@ociw.edu 

5 Southern Observatory for Astronomical Research, Casilla 603, La Ser- 
ena, Chile; sheathcote@noao.edu 

6 Department of Physics, Hofstra University, Hempstead, NY 11549; 
Stephen.Lawrence@hofstra.edu 



lines and a lack of early escape X and 7- rays) also exists. Fo r 
a complete accounting, see the review bv lArnett et alJ d!989h . 

Panagiaetal. (1987) reported narrow UV emission lines 
in International Ultraviolet Explorer spectra taken on day 
80, which they interpreted as arising from pre-existing CS 
materi al. An analysis of these spectra by iFransson et alJ 
( 1989) revealed a large overabundance of nitrogen, suggest- 
ing both significant CNO processing by the progenitor, and 
that a reasonable fraction of the progenitor's hydrogen enve- 
lope had to have been ejected to reveal such material. This 
argues that the star first evolved into a RS G before return- 
ing to the blue in a so-called "blue loop" (f Hofmeist er et alJ 
[1961. Both RSG JHurnnhrevs & Davidsonll1978l) and BSC, 
( Humphreys & McElrov 1984) stars are known within the 
LMC, thus this evolutionary sequence is consistent with 
known populations. 

The first res olved images of this CS material were taken 
on days 75 (Crotts Kunkel & McCarthy! 119891) and 1027 
JWampleret al.lll990al) . showing a central, elliptical struc- 
ture, s urround ed by two outer loops llCrotts et al.l Il989t 
Wampleretal. 1990b). The geometry of this material was 
the subject of some debate, remaining poorly understood un- 
til Hu bble Space Telescope {HST) imaging ( Burr ows et 
119951 and references therein) clearly revealed a dense cir- 
cumstellar equatorial ring (ER) flanked by two larger outer 
rings (ORs), as shown in Figure ^ To explain how a 
mass-losing supergiant produced this nebula, most authors 
invoke the interacting stella r winds scenario (Kwokl ll982t 
iBalick. Preston & "ick e 1987), in which colliding winds from 
the BSG and RSG mass-loss stages conspire to create three 
ring-like overdensities within the CS environment (CSE). 
These were ionized by the SN light pulse and are observed 
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FIG. 1.— HST WFPC2 negative image of a 4"5 X 5"0 field surrounding 
SN 1987A, taken in F656N (Ha and N II). North is up, east is left. The ER 
(central ring) surrounds the ejecta from the SN (center), and is flanked by the 
north and south ORs (NOR and SOR, respectively). The companion stars 2 
and 3 i Walker & Suntzeff 1 990 ) are positionally coincident along the line of 
sight. To increase the display range, the ER has a separate color stretch. 



today through recombination cooling. 

Since the SN lightcurve was strongly peaked for a fi- 
nite duration, and the progenitor is surrounded by a struc- 
tured CSE, we expect that the light pulse will scatter off of 
surrounding reflective material (e.g., dust), producing "light 
echoes" that a re observable long after the SN ha s fade d 
JCouderclll939t). iBond et al l Sl9S% . iCouch & Mairnl fl989l) . 
and Crotts & Kunkel ( 1989) discovered a light echo at ^9" 
whose position agrees well with the predicted location of a 
contact discontinuity (CD) between the RSG outflows and 
the interstellar bub ble formed by the early-type blue MS 
progenit or's win ds (Chevalier & Emmering 1989 ). Crotts et 
al. (e .g., ICrotts & Kunkelll 19911 ICrotts. Kunkel & Heathcotd 
1 19951 hereafter CKH95) have continued to search for light 
echoes in a regular campaign using emission-suppressed con- 
tinuum imaging. They have found echoes indicative of the 
structure of the three-ring nebulosity, forming an equatorial 
waist and sides of an hourglass, yet ending abruptly near the 
ORs with no capping surface. Outside this inner nebula, they 
found four additi onal feature s : (1) a reconfirmation of the 
CD light echo of lBond etal] il98% at 9-15"; (2) a sheet 
along the equatorial plane of the ER, bisecting the nebula; 

(3) diffuse flux from the hourglass nebula to the CD; and 

(4) a discontinuity i n this nebula, called "Napolean's Hat" 
( Wa mpleret al.ll990al) . 

Observation and modeling of this circumstellar environ- 
ment provide a unique opportunity to reconstruct the mass- 
loss history of the progenitor. The SN 1987 A nebula is one ex- 
ample of a much larger class of bipolar outflows, which occur 
in nearly all evolved stars (e.g., luminous blue variables, B[e] 
stars, AGB mass loss, bipolar and asymmetric planetary neb- 
ulae). This "stellar paleontology" is of immense importance 
for understanding bipolar-outflow and bipolar-nebula forma- 
tion mechanisms, not just of this object or in intermediate-to- 
high mass stars, but in stars of all masses. 

We summarize in this paper an observational effort, at the 
highest sensitivity and resolution to date, to understand more 
fully the entire circumstellar environment of SN 1987 A, and 
to recreate a more complete history of the progenitor's mass- 
loss, as revealed by light echoes in 16 years of optical imag- 
ing. After introducing echoes in <j2] we summarize the data 
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FIG. 2. — Cartoon schematic of a scattered-light echo. Echoes appear as 
arcs or rings on the plane of the sky, with a one-to-one mapping between 
their 2-D image positions (right) and their 3-D locations in space (left). N.B. 
Distances are not to scale, since D> z- 



and reduction in |3] These data are used to build a complete 
model of all material within ~30 ly of the SN, the 3-D anal- 
yses of which are summarized in 5j4] followed by the dust 
properties in ^5j Full details of these sections are presented in 
Sugerman et al. (2005, hereafter Paper II ). 

With a much greater volume of the progenitor's CSE re- 
vealed, we find that no extant model adequately explains 
how that environment was formed. In |6j we review and 
critique formation scenarios for this system, including the 
progenitor's MS wind, interacting supergiant winds, wind- 
compressed disks, and the influence of a binary companion. 
We use the light-echo data to constrain the parameter space of 
these models, and we qualitatively sketch a hybrid formation 
sequence for the entire CSE in It is our hope that the re- 
sults presented here will serve as a significantly-improved set 
of constraints for hydrodynamic models of the CSE forma- 
tion, for stellar-evolution models of the progenitor, and for 
addressing the general question of asymmetric and bipolar 
stellar outflows. 

2. LIGHT ECHOES 

When a light pulse is scattered by dust into the line of sight, 
an observable echo is produced, provided the pulse is suffi- 
ciently luminous and the dust sufficiently dense. An echo 
observed a given time after the pulse must lie on the locus 
of points equidistant in total light travel from the source and 
observer, that is, an ellipsoid with known foci. This simple 
geometry, shown in Figure [2J directly yields the line-of-sight 
depth z (and hence the three-dimensional, or 3-D, position) 
of an echo, uncertain only b y the assumed d istance D to the 
source. For fl» r, one finds dCouderd 19391) 



p z ct 
2rt~ 2" 



(1) 



where p = r sin 9 is the distance of the echo from the source in 
the plane of the sky, 9 is the scattering angle, and t is the delay 
between observing the echo and unscattered light pulse. 

A complete discussion of light echoes, including single- 
scattering models and the observability of echoes around a 
wide var iety of cataclysm ic and other variable stars, can be 
found in Sugerman (2003). The relevant light echo model, 
dust and outburst properties for SN 1987 A are summarized 
in §?? of Paper II. In general, echo brightness scales directly 
with dust density, inversely with distance, and decreases non- 
linearly with increasing scattering angle. 

3. DATA AND REDUCTION 

3.1. Observations 

SN 1987 A has been observed many times per year since 
day 375 (after core collapse) in a ground-based campaign 
to monitor the appearance and evolution of its light pulse 
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as it illuminates circumstellar and interstellar material. To 
minimize confusion from sources of nebular-line emission, 
four specially-selected continuum bands were used for the 
monitoring campaign, centered at 4700, 6120, 6880, and 
8090A. So me data from this campaign have previously bee n 
reported by ICrotts et alJ Jl989h : ICrotts & Heathcotd d!991h : 
ICrotts &Kunkell (119911) : a nd CKH95 in stu dying the circum- 
stellar environment, and in IXii. Crotts & Kunkellj1994[T995h 
and |Xu & Crotts (1999) for interstellar echoes. Seventeen 
epochs of imaging with high signal-to-noise, good seeing 
(arcsec or better), and sufficient resolution to resolve most 
crowded stars exist between 1998 Dec and 1996 Jan. These 
data are supplemented with archival HST Wide Field and 
Planetary Camera 2 (WFPC2) imaging taken between 1994 
and 2001. Full details of the observations are given in §?? of 
Paper II. 

3.2. Data Reduction 

The complete data-reduction is presented in §?? of Pa- 
per II. To search for light echo signal, we employ a cus- 
tomized implementatio n of th e PSF-matched difference imag- 
ing (Tomanev & Crotts 1996) IRAF package difimphot, fully 
described in §?? and Appendix A of Paper II. A subset of the 
PSF-matched difference images are shown in Figure |3j see 
§?? of Paper II for additional data and discussion. The in- 
nermost echoes were particularly bright up to day 1469, and 
are washed out at the brightness-stretch needed to display the 
outer echoes. In such cases, a 12" inset centered on the SN 
is shown at the top right of each panel, scaled in brightness to 
better-resolve the innermost echo signal. 

3.3. Echo Measurement and Visualization 

Echoes were detected and measured in difference images 
by fitting a series of convolved moffats to radial surface- 
brightness profiles of arclength 10° centered on the SN. See 
§?? of Paper II for complete details. The centers of all moffats 
have been marked in the difference images (Fig.|3jl, and color 
coded to compare echoes believed to be physically associated. 

During their detection, echoes were categorized as belong- 
ing to one of the three expected structures: (1) a circumstellar 
hourglass-shaped nebula reported in CKH95, (2) Napoleon's 
Hat (Wampler et al. 1990a), or (3) th e discontinuity between 
the pr ogenitor's RSG and MS winds ( C hevalier & Emmering l 
1989). In Figure Rl these have been color coded red, green, 
and blue, respectively. However we also found additional 
echoes between the hourglass and Napoleon's Hat, arbitrarily 
colored cyan for those looping to the north of the SN, or gold 
to the south. In later epochs, many of these southern echoes 
have positions that are more consistent with a counterpart to 
Napoleon's Hat, and have been marked purple. 

That diffuse structure exists between the inner hour- 
glass and Napol eon's Hat has been previously reported by 
ICrotts & Kunkell d!991h and CKH95. These are consistent 
with the cyan and yellow echoes identified in the same re- 
gion of the current data, which now map out a much more ex- 
tensive volume of the innermost CSE. Contact-discontinuity 
echoes have been previously reported in observations only up 
to day 1028, while we are able to trace this structure through 
day 3270. Napoleon's Hat was imaged as ear ly as day 850 and 
up to day 1650 by Wang & Wampler ( 1992), while we detect 
its signal from day 659 to 2095. We believe this is the first 
discussion of the southern-counterpart (purple) to Napoleon's 
Hat, which was illuminated between days 1787-3270. 



3.4. Echo Visualization 

To study echo positions in 3-D, we have custom-written 
graphics software to perform simple renderings, in which the 
positions of the echoes are transformed to allow viewing from 
any angle. The method is fully explained in §?? of Paper II, 
and briefly annotated in the caption of Figure^ which shows 
an illustrative example. 

The region of space probed by a light echo is an unusual ge- 
ometric function, and we are generally unfamiliar with view- 
ing structures by their intersection with parabloids. To fa- 
cilitate the understanding of this mapping, we have rendered 
echoes from a bipolar nebula in Figure |4] These are intended 
to guide the reader in visualizing such intersections, and in 
translating from 2-D data on the sky (left column) to 3-D po- 
sitions in real space (right columns). This also shows that the 
geometry of a structure can only be deduced from multiple 
echoes. See §?? of Paper II for more detailed discussion, and 
a map of the parabolae corresponding to all epochs of data. 

Echo positions in 3-D are computed from their epoch and 
2-D image positions using equation Q, assuming the echo 
center is produced by the light pulse maximum at day 87. 
We adopt a distance to the SN of 50 kpc, a com mon average 
between many of the derived LMC dista nces l lGould & Uzal 
1998: lFeasll999URomaniello et all20"00l 

The 3-D positions of all echoes are rendered in Figure [5] 
The top, middle, and bottom rows show the contact disconti- 
nuity (CD), Napoleon's Hat (NH), and the inner circumstellar 
(CS) material, respectively. Colors correspond to those in- 
dicated in Figure [3] except for the yellow points in the CD, 
which were identified in WFPC2 images very close to a bright 
star cluster northwest of the SN (PA. 300°, p = 15"). We dis- 
tinguish these data because large residuals from the closely- 
spaced, bright stars made it impossible to verify the echoes in 
ground-based images. 

4. A COMPLETE PICTURE OF THE CIRCUMSTELLAR 
ENVIRONMENT 

A rigorous analysis of the 3-D structure of the light echo 
data is presented in Paper II, §??-??, however we describe the 
results in the following subsections. We begin with the rea- 
sonable assumption that the echoes lie on continuous surfaces, 
only subsets of which have been illuminated. Echo brightness 
decreases with increasing scattering angle and distance from 
the illuminating source, and increases with dust number and 
column density. That echoes were not observed in some re- 
gions between the parabolae in Figure|5]need not imply lack 
of scattering dust, but may result from unfavorable geometry, 
or limitations in the data quality. 

4. 1 . Contact Discontinuity Echoes 

When examined in spherical-polar coordinates, the CD 
echoes (Fig. [Sja— c) appear to trace two different geometric 
shapes. Points east and west of the SN form a shell-like fea- 
ture at roughly constant radius, while points south of the SN 
(P.A. 150°-210°) form a radial feature at roughly constant in- 
clination. In both cases, the points appear to lie on a surface 
with some degree of symmetry, such as a spheroidal or coni- 
cal structure. Since only a partial cross section of the CD neb- 
ula has been illuminated, finding its inclination through nu- 
meric minimization of a merit function (e.g. the best-fit cone 
or line; see Appendix B of Paper II) was not highly successful. 
We supplemented those results by testing how the echo points 
compared to geometric figures observed at different orienta- 
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FIG. 3. — Selected 50" X 50" difference images. North is up, east is left, and major ticks mark 10". The position of the SN, Stars 2 and 3 are marked with 
yellow stars. The inset shows the central 12" at a different color stretch to resolve the innermost echoes. Echoes are marked by colored circles, explained in SI3.3I 
(a) 6067 A image from day 750. (b) 6120A image from day 1028. (c) 6120A image from day 1469. (d) 612A image from day 2095. (e) WFPC2 image from day 
2769. (/) 612A image from day 2874. 
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FIG. 4. — Examples of 3-D renderings of light echoes in a simple bipolar 
nebula, depicted at top, which is inclined 45° to the line of sight (to the right). 
Each lobe is 5 ly in radius, with the centers separated by 8 ly. Five echo 
parabolae are drawn, each occuring at the time (in years) indicated next to 
each curve. Below, renderings show the light echoes for the first two parabo- 
lae (first row), third parabola (second row), last two parabolae (third row), 
and all parabolae (last row). Axes are defined as z toward the observer and 
x - y in the plane of the sky, with x increasing west and y north. Axis la- 
bels indicate the positive direction. The origin is set at the SN, and is always 
marked by a black dot. Major ticks mark 2 ly, and the origin is indicated by 
the longest tick along each axis. The coordinates and axes have been given a 
slight perspective transformation (Appendix C of Paper II). Points are shaded 
using simple ray tracing, and larger points are closer to the observer. Left 
Column: Face-on view (the plane of the sky). Middle Column: Oblique view 
rotated 45°. Right Column: Side view from far to the east. Note that the 
geometry of the nebula, as revealed in light echoes, is only clear when many 
echoes well-separated in time are considered at once. 



tions, and conclude that the CD echoes lie roughly on a sur- 
face of rotation about a single axis inclined 40° south and 8° 
east of the line of sight. This axis may also be offset east of 
the SN by < 1 ly. The cross section gives some indication of 
an ellipticity of b/a = 0.95, with the major axis aligned toward 
north, however this is not a robust measurement. 
The data are transformed into a "primed" cylindrical co- 



ordinate system about this axis, designated by z' along the 
axis and tt' the radius. Figure [6] shows the resulting radial 
profile ir'{z')- Assuming these data trace a structure that (1) 
has rotational symmetry about z', and (2) is symmetric about 
the equator (z' = 0), an average radial profile of (it 1 ) versus z' 
can be constructed by binning the points along the axis. This 
(and standard deviations in each bin) is shown over the in- 
dividual data. Some circles located at tt' ~ 15 and z=4- 8 ly 
appear disconnected from the rest of the shell-like material. 
These points, which we call the "spurs," were treated sepa- 
rately when making the average radial profile. 

The probable CD structure is visualized by revolving this 
function about the symmetry axis, reflecting it about the equa- 
tor, and reinclining it to the favored orientation, as shown 
in Figure \l2h-b. Colors are the same as Figure |5j only the 
southern echoes have been colored red and green to distin- 
guish points marked with "x" and "+" in Fig. [6] Figure[T2T?-ij 
shows the complete structure in monotone grey, overlaid with 
the actual echo points from Figure [5] to show exactly which 
parts of the complete structure were sampled. 

The complete structure is fairly complicated, and can not be 
described by a single geometric function. Figure ITTb shows 
a scaled cartoon of the salient features. A "shell" lies along 
a prolate spheroid with a polar axis of 20 ly and equatorial 
axes of 1 1 ly. However, the ends of this spheroid have been 
drawn out into tapering, "radial" cones with opening angles of 
about 35°, extending from 16 to 28 ly from the SN. Embedded 
within this prolate structure are narrower, tapering cones (the 
"jet") extending from 10-26 ly from the SN, with an open- 
ing angle of about 20°, and a maximum radius of 3 ly. The 
"spurs" lie along a cylindrical annulus that smoothly encircles 
the CD. 

It is unclear whether the prolate shell is continuous in z' 
along its equator (dashed equatorial lines in Fig. II lit), since 
no echoes were observed from that region. Whether this is 
due to shadowing from material closer to the SN is addressed 
in §?? of Paper II. The northwestern- WFPC2 echoes are posi- 
tioned at z' ~ 4 ly along the axis, and at roughly 8 ly in radius. 
This could suggest the prolate shell is pinched at its waist to a 
smaller radius of 8 ly. However given the very limited spatial 
sampling of these inner echoes, there is little evidence that 
they lie on a uniform structure. Similarly, very little signal 
was detected from the spurs which may suggest that they are 
isolated clumps, rather than a uniform feature. 

4.2. Napoleon's Hat 

Figures show the 3-D positions of the NH echoes. 
The northern (green) points forming the familiar bow- shape 
for which this structure was named ( Wam pler et al J f 1 990a ) 
appear to lie on a thin, cylindrical shell that is completed 
by the newly-discovered southern points. To avoid confu- 
sion between this and previous work, we will keep the name 
"Napoleon's Hat," by which we refer to the entire ensemble 
of echoes. Despite the temptation to call the southern points 
"Napoleon's Collar," we will simply refer to these as NH- 
north and NH-south, when needed. 

Two nearly-horizontal features protrude from this shell to- 
ward the observer, appearing horn-like when viewed from 
above (panel f). In the following discussion, we refer to these 
as the "horns," to distinguish them from the rest of the NH ma- 
terial. These are the echoes seen in early imaging that bridged 
between NH-north and the larger-radii CD echoes (Fig.|5|l. 

Viewed in spherical-polar coordinates, the NH echoes are 
consistent with a cylinder or hourglass. We measured the 
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FIG. 5. — Rendered views of all light echoes identified in Fig. [3] Points have been rendered using the method described in SI3.4I also see Fig. [3] The left column 
shows observed views in the plane of the sky, the middle column shows views from the side (far to the east) and the right column shows views from the top (far to 
the north). Parabolae from the earliest and latest epochs at which echoes were observed are indicated in the side and top views. Point colors correspond to those 
in Fig.|3]except for the yellow points in the top row, which denote echoes found in WFPC2 images which could not be resolved from the residuals of nearby 
bright stars in ground-based images. Note that the field-of-view changes between rows, but major tick marks always denote 2 ly. Top row: Contact discontinuity 
echoes. Middle row: Napoleon's Hat and southern counterpart. Bottom row: Inner circumstellar hourglass echoes. 



shape and inclination of such structures by fitting a biconical 
frustum to all points, as explained in Appendix B.4 of Paper 
II. The best fit is an hourglass with elliptical cross section, 
oriented 40° south and 7° east of the line of sight. The semi- 
major and minor axes are 4.6 and 3.8 ly, with opening angles 
at the axes of 28 and 40°, respectively, and the major axis is 
rotated 92-101° east of north. That the cross-section is ellip- 
tical is unexpected, and we reverified this by directly examin- 
ing the data looking down the axis of inclination. The best-fit 
ellipse (Appendix B.2 of Paper II) to the resulting distribution 
has negligible offset from the origin, and b/a = 0.82 ± 0.02 
with the major axis rotated 103° ± 3 east of north. 

Removing this ellipticity from the radial values of %' as a 
function of z' yields the distribution shown in Figure Also 



shown is the average radial profile of the NH echoes, gener- 
ated by binning the points along the inclination axis and com- 
puting the average radius about that value. As with the CD, 
we visualize the probable structure containing these echoes by 
reflecting this average radial profile about the equator and re- 
volving it around the axis. The result is an inclined hourglass, 
as rendered in Figurefl2fe-/i. 

An examination of Figure [T2h shows that very little of the 
hourglass was actually probed by the horns. We show in §?? 
of Paper II that such limited illumination was a natural result 
of the scattering geometry of this complete structure. The 
hourglass therefore flares to a radius (or semi-minor axis if 
the cross section remains elliptical) of ~ 7 ly about 7 ly from 
the SN along the axis of symmetry. 
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FIG. 6. — CD echoes, transformed into cylindrical coordinates with z' along 
the axis of symmetry (40° south, 8° east) and 7r' measured radially from it. 
Eastern and western points are marked by filled and open circles, respectively. 
Southern points are marked as "x" except for a subset at smaller it' marked 
with "+". A small subset of WFPC2 echoes to the northwest (yellow points 
in Fig.[3J are marked with asterisks. Overplotted is the average radial profile, 
measured as the average position (and standard deviation) of points binned 
along z! ■ 




\A (iy) 

FIG. 7. — As Fig.|5| but for NH data, with the ellipticity removed. Over- 
plotted is the average radial profile, measured in bins along \z'\. Single points 
have fixed errors. 



The complete cross section is shown in the cartoon sketch 
in Figure[Tlt>. Of note are the radial features in, and about 1.5 
ly above/below, the equatorial plane. CKH95 reported a small 
set of echoes within 3 ly of the SN that are coincident with the 
ER plane, which they interpreted as evidence for an extended 
circumstellar equatorial disk. These NH equatorial echoes ap- 
pear to extend this disk to a radius of 5 ly, with a thickness of 
> 0.75 ly. We note that this structure is inconsistent with NH 
being a parabolic b ow- shock north of the SN, as proposed by 
Wang & Wampler ( 1992) and lWang. Dvson. & Kahnl i 19931) . 
This is discussed in greater detail in §?? of Paper II. 

4.3. Inner Circumstellar Material 

The inner 3" surrounding SN 1987A have been previously 
studied by CKH95 using a subset of the data presented here, 
in which the authors reported a double-lobed CS hourglass 
with a distinct cylindrical symmetry axis. Although the red 
echo points in Figure appear to lie along such a surface, 
we cannot immediately explore this feature in our data since it 
is unclear which extended-flux echoes (cyan and gold points) 
are associated with an hourglass. We use the work of CKH95 
as a springboard to help disentangle our data. 

Fitting a biconical hourglass to the CKH95 data, we find the 
structure is inclined 45° south and 8° east of the line of sight. 
This compares favorably to the orientation they reported us- 





(a) (b) 
FIG. 8. — CS points believed to lie along the inner hourglass. 
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FIG . 9 . — As Fig.|^] but showing the CS data. Inner hourglass points (red in 
Fig.l5t-t) are marked with "x", extended-flux points (cyan and gold in Fig. 
Bl?— fi with circles, and extended flux points belonging to the hourglass with 
squares. The average radial profiles to each structure are overplotted. Note 
that there are many binned averages in this plot, which are discussed in the 
text. 



ing a minimization-of-scatter estimator. Adopting this as a 
preliminary orientation, we transformed our CS echo points 
into cylindrical coordinates about this axis, which revealed an 
additional subset consistent with those initially identified as 
part of the hourglass. The full set of probable hourglass points 
is rendered in Figure [8] As can be seen in panel (b), only a 
small arc along the southernmost limb of of the southern lobe 
was probed by the earliest observed echoes. 

A combination of fitting biconical hourglasses to these data, 
and ensuring that the inner waist of this hourglass reproduces 
the observed shape of the ER, yields the 3-D geometry of both 
structures. We find the hourglass is inclined 41° south and 8° 
east of the line of sight, with its axis shifted 0.1 ly west of 
the SN. The waist has a semi-major axis of 1.04 ly rotated 
9° north of east, with b/a = 0.94, and the frusta have half- 
opening an gles of 12. Note that t his is consistent with the 
findings by ISugerman et ail (|2002) that the nebula is offset 
slightly west of the SN. 

Deprojected by the above inclination, the ER is best-fit 
by an ellipse with major axis 0"82 rotated 9° north of east, 
bja = 0.98, and the centroid is shifted 19 mas west of the SN. 
This is the first direct measurement of the ER's orientation 
and deprojected geometry. Full details of these analyses are 
presented in §?? of Paper II. 

As in the previous two subsections, we incline the data 
about this axis, remove any ellipticity, and plot the cylindrical- 
radial profile in Figure [9] Extended flux points that are now 
associated with the hourglass are plotted as squares. A large 
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subset of the extended-flux echoes (circles) are concentrated 
within 1 ly of the equatorial plane (z' = 0), between 1 .5-4 ly 
from the SN, and appear to form a thick equatorial waist or 
belt circumscribing the hourglass. As in the previous two sec- 
tions, we measure the average radial profile of these points in 
bins along |z'|. There is much structure to these data, which 
required the fitting of six profiles. We again visualize the neb- 
ula by reflecting each profile about the equator, and revolving 
it around the cylindrical axis, as rendered in FigurefT^t-p. and 
sketched in cross section in FigurefTTl:. We identify five fea- 
tures in this data, discussed below. 

The innermost radial structure is the hourglass, colored red 
in Figure \\2\-p. This is fairly cylindrical (opening angle 
~12°), with a semi-minor axis of 1.2 ly, tapering close to the 
equator and opening just slightly at large z' ■ A very small 
number of points suggest that it also tapers around z' = 1 .5 -2 
ly and flares at z' = 0.9 and 1.4 ly, but these features are 
marginal. Consistent with CKH95, we do not see evidence 
for a "capping surface" to this structure, as would be expected 
from a peanut-like nebula. To distinguish this hourglass from 
the larger one containing the NH echoes, we refer to each by 
their association, i.e., the former is the "CS hourglass." 

In Figure [9] many of the extended-flux points (squares) lie 
just outside the average position of the hourglass. Since most 
of these points are part of the southern lobe, this could sug- 
gest it is wider than its northern counterpart. A more likely 
interpretation is that many of the southern-lobe points repre- 
sent dust just outside the densest part of the hourglass, which 
was itself illuminated prior to our earliest epoch. 

Immediately surrounding the CS hourglass is the thick 
waist of extended high-surface brightness, which we call the 
"belt" (colored green in Fig, \l2\-i). It extends ±1 ly along 
the z' axis, and from 1.5 to 2.5 ly in radius. Beyond the outer 
radius of this belt, echoes lie fairly well constrained to the 
equator, tracing a thinner (0.5 ly thickness in z') but extended 
equatorial plane (blue in Fig.ll2t-/") to an outer radius of n' = 4 
ly. The belt does not appear to taper smoothly into the thinner 
equatorial plane, but rather has a fairly sharp transition at its 
outer boundary. 

The rest of the CS material lies along one of two "walls." 
The inner wall is colored cyan in Figure fT2t-/'. extending 1.5 
to 2.6 ly along the axis, and may be considered to join with 
the hourglass around \z'\ = 1.5 ly. The outer wall, colored 
gold, bridges the belt and the inner wall at radii between 2.6 
and 3.4-4.0 ly. We note the conspicuous lack of echo signal 
between the hourglass and outer wall around \z? \ = 1 .4 ly, and 
question whether this is real, perhaps indicative of episodic 
mass loss, or simply a gap in our data. 

4.4. The Geometry of the Outer Rings 

While great attention has bee n plac ed on the geom- 
etry of the E R (Jako bsen etafl H99ll iPlait et alJ Il995l 
iBurrows et al.lll9 95l). rel atively lit tle work exists o n the ORs 
(see B urrows et al.ll995tlCrotts & Heafhcotd2000l) . Now that 
a 3-D map of the nebula exists, we can investigate how it can 
constrains the geometry of the ORs (full details are in in §?? 
of Paper II). 

The North OR (NOR) appears to lie along the outer edge 
of the north CS hourglass, while the South OR (SOR) seems 
to lie at the intersection between the southern hourglass and 
the belt, perhaps indicating an interaction between wind that 
created the hourglass and the pre-existing belt material. 

The 3-D positions of the ORs are rendered in Figure ITOl 
and appear reasonably planar. We approximate both rings as 
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FIG. 10. — Probable positions of the ORs (grey) and the planar, elliptical 
fits (black) to each ring, viewed (a) face-on, {b) 60° east, (c) 60° west, and 
(d) orthogonal to the inclination axis. 



TABLE 1 
Best-fit ellipses to the CS rings 



Ring 
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y'o 
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a 


b 






(ly west) 


(ly north) 


(ly) 


(ly) 


(ly) 


(°) 


ER 


0.015 


0.0 


0.0 


0.647 


0.98 


81.1 


NOR 


0.26 


0.04 


-1.36 


1.42 


0.94 


70.5 


SOR 


0.19 


0.06 


1.00 


1.59 


0.92 


-1.1 



planar ellipses, fit to each distibution of points viewed along 
the CS inclination axis. The fitted parameters are listed in 
Table [l] where (x' Q ,y' ) are the centroid offsets from the axis, 
Zq is the distance to the SN along the axis, (a,b) are the semi- 
major and minor axes, and <f> is the RA. of the major axis from 
north. For completeness, the parameters for the ER are also 
listed. These fits are also indicated in black in Figure^H The 
approximation to the NOR is quite good, but the SOR ellipse 
fails to intersect the points that pass just north and west of 
the ER. This is not suprising, since no single ellipse can fully 
reproduce the observed shape of the SOR unless it is non- 
planar. We have used these planar approximations in the echo 
renderings throughout this work. 

4.5. Summary 

The probable geometry and structure of the echoes within 
25" of SN 1987A are rendered in Figure^l ar >d summarized 
in Table [2] The final column gives the approximate total vol- 
ume of the structure containing each echo, found by integrat- 
ing the radial profiles in Figures [6] [7J an d HI along, and re- 
volved 2ir about, the respective inclination axes. 

The east/west and north/south inclinations are quite consis- 
tent, suggesting all circumstellar material shares a common 
inclination near i x = 40° south and i y = 8° west. Furthermore, 
the structures appear to join to form a single nebula, as fol- 
lows. The CS gas is nested neatly within the NH material, 
as shown in Figure fT2T. with the CS belt, equatorial plane, 
outer wall, and the waist of NH consistent with a single, uni- 
form thick waist, extending inward to the CS hourglass (Fig. 
112b). The outer edges of the NH hourglass almost reach the 
inner spur of the CD shell (Fig. 112b ). and indeed the echoes 
as viewed on the sky do connect at early times (Fig. [5ji. As 
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Equatorial 
^ plane 




^ (b) 

FIG. 11. — Cartoon sketches of the salient structures traced out by (a) the CD model from Fig.|6| (b) the NH model from Fig.Q and (c) the CS model from Fig. 
|9| Figures are to scale, with the length of the two orientation arrows equal to (a) 10 ly, (b) 4 ly, and (c) 2 ly. Solid lines trace the radial profiles, and the widths of 
the greyscale regions trace the scatter of points about the averages. 



shown in Figure [T2lw. no material was illuminated along the 
equatorial plane exterior to NH. We do not believe this re- 
gion was shadowed from the SN by material at smaller radii 
(^5]of this paper, or §?? of Paper II), thus it is unlikely there 
is any higher density gas outside the waist of NH. As such, we 
conclude the NH hourglass is the pinched waist of the CD. 

Figure ^] shows all circumstellar material, simplified ac- 
cording to the above arguments. The outer structure is peanut- 
like, extended along the poles to r < 28 ly, and narrowly- 
pinched at the waist at r ~ 5 ly. We will refer to this outer 
shell (traced by the CD and outer NH echoes) as the "Peanut." 
Roughly 35° from the equatorial plane are the spurs, extend- 
ing 3 - 4 ly out from each lobe, perpendicular to the Peanut's 
axis. Unlike the northwestern WFPC2 echoes, these spurs do 
appear at more than one location, but they are limited enough 
in extent that it is unclear whether they are uniform features 
that encircle the entire structure, or isolated clumps. 

Inside the equatorial region of the Peanut is a thick (~ 3.5 ly 
in radial thickness, 2 ly in axial length) annular ring (or belt), 
which extends inward until it terminates at the inner CS hour- 
glass. The extent of this waist in z' is unclear. Echoes were 
detected outside the CS hourglass at most equatorial (ir 1 ) radii 
between 1.6 and 4 ly, with a vertical (z') distribution extend- 
ing at least as far as the CS hourglass itself. Still, the bulk of 
material is within 1 ly of the equator. 

5. DENSITY AND MASS OF THE CSE 

To constrain the gas density and dust composition of the 
CSE, we apply a dust-scattering model (§?? of Paper II) using 
the surface brightnesses of all echoes discussed in the previ- 
ous section. The reader is referred to §?? of Paper II for the 
complete analysis, as only a brief summary is presented here. 

The CD echoes are best fit with large (a = 3.5A-2.0/im) 
grains consistent with an LMC-abundance of silicate and car- 
bonaceous dust. Outer NH data (r > 5 ly) are more consistent 
with Si-dominated large grains, but perhaps with a slightly- 
smaller maximum grain size than for the CD. The inner NH 
points, along with the CS echoes, are well fit by Si-only dust 
with a max = 0.2/im. The general trend is that inner echoes are 
better reproduced by smaller, Si-dominated dust, and with in- 
creasing distan ce, the grain sizes and C-content increases. 

Fischera. Tuffs. & Volk ( 2002) have modeled mid-IR emis- 
sion from collisionally-excited dust grains heated in the 
shocks between the CS gas and the forward blast of the 
SN. They find the emission is best explained by small (a < 



0.25/^m) grains with a Si-Fe or Si-C composition. The dust 
abundance is quite low, which they attribute to evaporation 
from the UV flash and sputtering in the shocked gas. This 
further constrains grain sizes to a < 0.25/im and excludes a 
pure-carbon composition. The dust in both this shocked re- 
gion and the CS hourglass was formed from material expelled 
very late in the progenitor's life, thus we expect the dust prop- 
erties to be similar. Indeed, our pre-SN CS dust model also 
favors silicate-dominated dust with grain sizes a < 0.2/im. 

That dust which formed at later times has a smaller carbona- 
ceous content can result from a change in surface CNO abun- 
dances over the star's late-stages of evolution, since carbon- 
rich envelopes create carbon-rich dust, while oxygen-rich en- 
velopes create silicate-rich dust. Such CNO processing is also 
inferred from early IUE spectra (Fransson et al. 1989) of the 
ER, which show nitrogen and oxygen to be overabundant with 
respect to carbon. 

The average dust density is plotted in Figure ^] Figure [T3l 
has been greyscale shaded to reflect the differences in density 
among the many echoing structures, with darker grey indicat- 
ing higher-density material. 

The innermost CS hourglass material has a relatively- 
constant gas density of «h = 2-3 cm" 3 up to r ~ 1.6 ly. Be- 
yond this position, the material splits into two distributions, 
where the steeper profile (r -4 5 ) traces the denser waist of the 
hourglass. Beyond r = 2 ly, the shallower (r~ 3 ) profile traces 
the belt and equatorial plane material, making a smooth tran- 
sition to the inner NH points around r > 3 ly. That the CS 
density varies more steeply than the inverse-square expected 
for a freely-expanding wind suggests the mass-loss mechani- 
cal luminosity {Mv^ xp ) increased with time toward the end of 
theRSG. 

Excluding a density enhancement between r = 4 - 5 ly, 
there is no evidence of a structural distinction between the 
CS and NH gas, which justifies the simplified model in Fig- 
ure Returning to the aforementioned enhancement, this 
higher-density material is located along the waist of the NH 
hourglass, marking an outer edge to the equatorial overden- 
sity. The constituents of the Peanut, the CD and outer NH, 
have constant density with radius, bounded by 0.02 < «h < 
0.2 cm" 3 , suggesting this structure marks a constant-density 
b oundary. 

iBurrows et alJ (11995ft argued that ;/ a CS hourglass exists, 
the fact that it is not observed in recombination (like the three 
rings) implies its density must be < 5 cm -3 . This density is 
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FIG. 12. — Rendered views of the probable structures containing the observed echoes. Top Row: The CD nebula. Second Row. The NH nebula. Third Row: 
The CS nebula. Fourth Row: The CS hourglass. For these top four rows — Left column: the face-on view of the complete structure; Second column: the western 
half, viewed from the east, showing a clear view of the interior; Right columns: face-on and side views of the complete structures in monotone grey, overlaid 
with actual echo points from Figure|5] For the bottom row — (q) The western halves of the CD and NH, and (r) NH and CS nebulae, viewed from the east, (s) 
The northern halves of the CD and NH, and (f) the NH and CS nebulae, viewed along the inclination axes from the south. See text for colors. 



consistent with our findings, and explains why the hourglass 
exists but is optically invisible. 

In §?? of Paper II, we present an unbiased measurement 
of density by searching for echoes at all positions within each 
difference image. These reveal the region between the CS and 
CD structures to be filled with diffuse material («h ^ 0.03), 
while little structure is seen outside the CD. According to 
Figure^] the CD and NH form a peanut with a narrowly- 



pinched waist. If the equatorial dust density is high enough to 
cause significant extinction, material beyond NH will be ef- 
fectively shadowed from the SN light pulse, in which case no 
echoes would be seen. In §?? of Paper II, we show that only 
a fraction of a magnitude of extinction can be expected from 
known equatorial material. As this is insufficient to cause pro- 
nounced shadowing, the narrow waist appears to be a genuine 
boundary. 
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TABLE 2 
Summary of Echo Geometries 



Struc- 
ture 


Inclination" 
South East 

(°) (°) 


Roll 

o 


Cross Section 13 
b/a A.v Ay 

(ly) (ly) 


Shape 


r 

(ly) 


z' 
(ly) 


Geometry' 

tt' 

(ly) 


Att' 

(ly) 


V M 
(ly 3 ) 


CS 


41° 


-8° 


-9° 


0.94 


<o.i 


0. 


Hourglass 


1-2.9 


0.3-2.5 


0.8-1.6 


0.5-1 


20 
















Belt 


1.5-2.8 


0.-1.0 


1.5-2.6 


1.5-2.0 


40 
















Walls 


1.8-4.2 


0.8-2.7 


1.8-3.6 


2.0-3.0 


120 


NH 


40° 


-7° 


12° 


0.82 


0. 


0. 


Hourglass 


3-11 


0.7-7.5 


2.8-7.5 


2-4 


1300 


CD 


40° 


-8° 




>0.95 


<-l 


<-0.1 


Prolate 


10-28 


4-28 


0-16 


4-5 


1.2 x 10 4 



""Roll" is the counter-clockwise rotation about the z axis, also the RA. of the minor axis after removing the southern inclination and eastern 
rotation. 

b Axq and Ayg are the western and northern offsets of the symmetry axis from the SN, which is at x = 0, y = 0. 

c r is the spherical-polar distance from the SN. 7r' is the cylindrical-polar radius measured from the inclined axis z', A-7r' is the approximate 
width along 7r'. V, ot is the inferred total volume of the structure. 




FIG. 13. — Revised ca rtoon showing the simplified CSE suggested by com- 
bining the panels in Fig. II II Structures are shaded to indicate density, which 
increases with greyscale. Figure is to scale, and orientation arrows indicate 
10 ly. 



The total mass of gas and dust within the CSE is estimated 
from the average density and volume of each structure, as 
given in Figure [H] and Table |5] The volume interior to the 
CD is estimated a few ways: as the volume of the ellipsoid 
enclosed by the CD, as the volume of two cones enclosed by 
the quasi-linear profile, and as the integrated volume under 
the same profile (Fig.|6}. All three methods give a consistent 
result of about 6500 ly 3 . The calculated masses are listed in 
Table |3 yielding a total mass in all structures of ~ UMq 
(assuming an LMC gas-to-dust mass ratio of 400-600). It is 
very difficult to quantify all the sources of error, but we be- 
lieve the tendency will be to underpredict the masses of the 
structures, given conservative measurements of their volume 
and the unknown extent to which dust was not observable, 
due to unfavorable geometry, poor data quality, confusion, or 
possible shadowing. 

6. PREVIOUS FORMATION MODELS 

As discussed briefly in §l\ the progenitor's pre-explosion 
spectral type, as well as the early SN evolution, argue Sk 
-69° 202 exploded as a BSC Evidence of pre-existing CS 
material, and nitrogen enrichment within that material, imply 
the star first passed through a RSG phase, executing a "blue 
loop" to return to a BSG prior to core collapse. 

Mass loss can explain this blue loop, however only 
if the star loses either a small fraction or nearly all 
of its hydrogen envelope (IWooslev. Pinto &Ensmai] fl988: 



Wooslev, Pinto & Weaveri ll988l iNomoto et alM988l). While 
the latter option is ruled out by observations (Dopita! ri988l 
Walb orn et al . 1989), even low mass loss can qualitatively ac- 
count for th e observed CSE via interacting stell ar winds (or 
ISW lKwokl[l98l iBalick. Preston & Ickdll987l) . Originally 
proposed to explain planetary nebulae, a slow, dense wind 
from the evolved star (or RSG) is overtaken by a fast, ten- 
uous one, accelerated from the hot, blue core (or BSG). In 
particular, inhomogeneities in the winds, such as equatorial 
overdensities, shock focus the winds into characteristic bipo- 
lar shapes. 

Unfortunately, the formation mechani s ms of such over- 
densities are poorly understood. iMorrisI Jl98ll fl987i) pro- 
posed that a binary companion can create an overdensity 
in the primary wind within the pair's orbital plane, while 
iSoker & Liviol dl99l argue that this configuration can drive 
high-velocity jets out of that plane. Two models exist for 
the creatio n of a disk by single stars . A rotation-induced bi- 
stability (Lamers & Pauldrach 1991) produces a larger mass- 
loss or lower wind-velocity at the equator than the poles, re- 
sulting from the high sensitivi ty of line-d riven winds to opti- 
cal depth. Biorkman & Cassinelli ( 1993) have introduced the 
model of a wind-compressed disk (WCD), in which the orbital 
streamlines of gas launched from the northern and southern 
hemispheres of a star collide, producing an equatorial over- 
density. 

Many of these mechanisms have been invoked to explain 
the formation of the three-ring nebula, which has been the 
most prominent feature of the CSE. In the following subsec- 
tions, we review these models, and compare their predictions 
to the complete nebula revealed by light echoes. In short, we 
find that no extant model can adequately reproduce this sys- 
tem, however a number of features among them may offer 
some insight into the progentor's probable mass-loss history. 

6.1. Interacting Stellar Winds Models 

iLuo&McCravl ^9 lb) first invoked the ISW model to 
explain the emerging picture of an hourglass-shaped nebula 
contai ning the three rings. Us ing the thin-shell approxima- 
tion JMac Low & McCravl ll988l) in which all swept-up ma- 
terial lies in a thin shell with an isobaric interior (effectively 
a momentum-conserving snowplow), they simulated the col- 
lision of a fast, tenuous BSG wind with an older, slow and 
dense RSG outflow. Since, this colli sion has been s tudied 
with 2-D h ydrodynamic sim ulations bv lBlondin & Lun davist 
( fl99l and lMartin & Arnet^(U995l) . 
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FI G. 1 4. — Average gas densities that match observed echo fluxes to the dust-scattering model. Point colors correspond to those in the two right columns of 
Fig. 1121 CS material is marked by circles, NH by squares, CD by triangles, and smaller symbols denote measurements with S/N < 2. Dotted lines are the best-fits 
through the upper and lower CS loci, corresponding to n oc 1 (upper) and n oc r (lower). Horizontal dashed lines delineate the rough density boundaries of 
the CD gas. 



A number of results are similar between all of these mod- 
els. (1) A bipolar or hourglass nebula can be formed via this 
ISW framework that roughly matches the size and expansion 
velocities of the rings. (2) To match these constraints, an ex- 
traordinary asymmetry is required in the RSG mass loss, pro- 
ducing a significant equatorial density enhancement over that 
at the poles. We note that this conclusion is a result of the 
simulations, not an input condition. (3) The wind parame- 
ters necessary to produce this nebula are consistent with av- 
erage values of other stars, with the RSG characterized by 
a mass loss Mr = 6 x 10~ 6 - 10~ 5 Mq yr -1 , wind velocity 
v R = 8- 10 kms" 1 , and lifetime r fi < 10 6 yr, and the BSG 
byM B = 1.5 x 10 _7 -3 x 1O" 6 M yr" 1 , v B = 300-550 km s" 1 , 
and Tg = 10 4 yr. 

As predicted above, we do find an extended equatorial over- 
density in Figure ^] However, compared to this newly- 
revealed CSE, the models all have significant shortcomings. 
Foremost, nearly all the predicted CS nebulae are quite flat- 
tened and oblate (mushroom shaped), inconsistent with the 
prolate, hourglass shape we have mapped out. To generate 
roughly the right shape, the equatorial expansion velocities 
become far too large, while matching these velocities can lead 
to RSG mass-loss rates that are too high. Although these mod- 
els can explain the formation of the ER, they fail to explain 
the existence of the ORs. This may be linked to the fact that 
almost all model hourglass nebulae have solid capping sur- 
faces, which we do not observe. Finally, the adoption of a pre- 
existing RSG asymmetry is entirely ad hoc, and these models 
make no attempt to explain how such an asymmetry formed, 
nor to model what geometry its interface with the MS bubble 
should be. 

6.2. A Wind-Compressed Disk Model 

ICollins et alJ (Q999) have attempted to address one short- 
coming in previous work by assuming that the RSG wind 
forms a WCD from which the density asymmetry is derived. 
In contrast to the above work, they also allow for a veloc- 
ity asymmetry in the RSG wind, as defined by the disk. A 
WCD is characterized by the rotation parameter = V rot /V cr it, 
the ratio of the stellar rotation rate to the break-up speed. 
For rotations less than a threshold value £1 < the star 
forms a wind-compressed zone with an equator-to-pole den- 
sity contrast of 3-10, while for il > f2,/, a WCD forms with 
density contrasts ~ 100. By varying the WCD parameters, 
ICollins et all {l999) varied the RSG-wind environment sur- 
rounding the star, establishing the initial conditions for a BSG 



wind. This they addressed usin g self-similar solutions of 
wind-b lown bubbles from Dwark adas. Chevalier. & Blondinl 
(1996). 

For their fiducial model (a RSG with 20M Q , which is too 
high from stellar-evolution arguments), they find ~ 0.3, 
which produces a density contrast of 25. A variation of the 
rotation parameter between = 0.3-0.36 produces peanut- 
shaped nebulae with equatorial expansion velocities of 8-10 
km s -1 and a maximum (cylindrical) lobe radius roughly twice 
that in the equator (consistent with the diameters of the ORs 
and ER), for RSG wind parameters Mr ~ 10" 5 M yr" 1 and 
Vr = 20 km s" 1 , and BSG parameters M B ~ 10" 7 M Q yr" 1 and 
vb = 400 km s" 1 . A WCD is not required to produce a peanut, 
as a dense equatorial zone is sufficient. Again, the peanut 
shape does not fit the observed CS hourglas s, however they 
addre ss this by noting that in simulations ( Frank & Mellema 
the BSG wind can be shock-focused into the polar 
direction, resulting in a substantially-more prolate structure 
than that predicted from a self-similar solution. 

A few results stand out from this work. (1) The WCD is a 
straightforward explanation for the observed equatorial den- 
sity enhancement. (2) The authors do not find the need to 
invoke unusually low BSG wind momentum or RSG velocity, 
although the RSG mass-loss rate is much higher than we de- 
duce in 37.21 (3) A single RSG star cannot rotate at O > 0.25 
withou t the influence of a binary companion. Colli ns et alJ 
( 1999) estimate that a lower limit of 0.6M Q for the compan- 
ion is sufficient to produce the necessary rotation speeds for 
this model. If rotation is to be invoked as a formation mecha- 
nism, then the progenitor had to have been a binary system. 

6.3. Binary Models 

B inarity has been invoked bv | Posd iadlowsk i~& Jossl l[l989) 
and Podsiadlowski. Joss. & Ranpaport ( 1990) among others 
to explain the observed characteristics of the SN. A com- 
panion could have been engulfed by Sk -69° 202 during a 
common-envelope phase; however, for lower mass primaries, 
the secondary can cause d ynamical instabilities that eject the 
primary envelope ( Rasio & LiviolfT996l) . a process that is not 
s uggeste d from the SN evolution. 

ISoiceil £1999) invokes the influence of a binary compan- 
ion to produce the observed three-ring nebula, but his results 
are inconsistent with the ensemble of d ata presented here. 
Podsiadlowski. Fabi an. & Stevens! i 19911) suggest that the col- 
lision of winds between the primary and secondary spread out 
the wind into a truncated double cone. This is a novel forma- 
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tion mechanism for the hourglass feature we observe; how- 
ever, their geometric model further predicts that Napoleon's 
Hat lies on the northern-half of the cone, which we now know 

to be false. 

iLlovd. O'Brien. & Kahnl i 19951) suggest a variant on the 
double-cone scenario, in which the BSG wind sweeps up the 
base of the cones, leaving the end surfaces as the ORs. Fur- 
thermore, dynamical friction causes the secondary to coalesce 
into the primary, ejecting the envelope, thereby forming the 
ER. Again, this model is attractive as it suggests that a sin- 
gle mechanism formed all three rings. A testable prediction is 
that extensive RSG double-cone material should exist at and 
beyond the ring positions, however we see no evidence for 
such material in light echoes. If binarity is to explain the for- 
mation of the three rings, these models must be revisited to be 
consistent with the new picture of the complete CSE. 

6.4. Ionization 

All the above models have neglected the effects of ioniza- 
tion on the evolution of ISWs. The progenitor MS star is ex- 
pected to emit an ionizing flux of S = 2 x 10 48 7 s" 1 , wh ile the 
B3 I BSG star should emit S = 4 x 10 45 7 s" 1 dPanagial 19731) . 
both of which are sufficient to drive an ionization front ahead 
o f the outer shock. 

IChevalier & D warkadas (119951) have addressed the role of 
ionization in the to explain the reappe arance of radio emis- 
sion three years after outburst (Manchester et al. 2002). Their 
model of a pre-SN H II region suggests that the BSG had 
a mass loss of only M B = 8 x 10" 8 M Q yr" 1 with v B = 450 
kms" 1 . Furthermore, the ionization front that formed this 
H II region would move rapidly through the lower-density po- 
lar regions of the hourglass, and could have broken out of 
the lobes, thereby eliminating any capping-surface to the CS 
hourglass. This is a particularly interesting scenario for ex- 
plaining the observed geometry of the CS hourglass, as well 
a s the evolution of the SN remnant across the spectrum. 

Mever ( 1997) has proposed that an ionization front in the 
BSG phase induces hydrodynamic motions in the pre-existing 
asymmetric RSG wind. These motions create a latitudinally- 
dependent density profile that peaks about 50° abo ve the 
equato r, w hich he identifies as the out er rings, however lSokerl 
(1999) and Crotts & Heathcotel (120001) have offered a few ob- 
jections to this scenario. Fully-radiative hydrodynamic simu- 
lations including ionization are necessary to test these ideas. 

7. CONSTRAINTS ON THE EVOLUTION OF SK -69° 202 

We now investigate what constraints the geometry and den- 
sity of the CSE can place on the progenitor's evolution. 

7.1. The Main Sequence Bubble 

We begin by discussing the possible mass-loss scenarios for 
the MS progenitor, since the impact of these winds on the 
surrounding ISM for ms the initial condition s for subsequent 
mass-loss evolution. IChevalier & Emmering l (119891) note that 
the wind should have blown a bubble > 10 pc in radius, 
whose interior homogenized as the star evolved onto the RSG 
branch. 

Let us assume th at Sk -69° 202 was a zero-age MS type 09 
V with 18-20M Q ( IChevalier & EmmeririglU98 9l luminosity 
10 4 9 L Q and radius 8/?© ( Pana giall973l) . yielding a MS mass- 
loss rate of Mms ~ 10~ 7 M Q yr" 1 ( Nieuwenhuiizen & de Jager 
1990). The interaction of a lo ng-lived stellar wind with 
the ISM was first investigated bv ICastor. McCrav. & Weaver] 



TABLE 3 
Mass of the CSM 



Structure 


M/Mq 


CS Hourglass 


0.04 


CS Belt 


0.04 


CS Walls 


0.06 


NH Waist 


0.06 


NH Walls 


0.07 


Intra-CD 


0.2 


CD 


1.2 


Total 


1.7 



i 19751) and lWeaver et alJ ill 9771). with refinements p r opose d 
by, among others. iMcKee. van Buren. & Lazarefn (1984), 
Koo & McKee (T992a.bl and iD'Ercold (119921) . The 
wind is characterized by its mechanical luminosity L w = 
jMv^, where the terminal velocity has an upper limit 
around < 1500 km s" 1 for the assumed progenito r 
properties (Kudritzk i et alJI 19891 Earners & Cassinellilll999T) . 
These yield L36 = 0.07. wher e L36 = L w /10 36 ergs s" 1 . 
Wooslev. Pinto & Weaver ( 1988) suggest the mass loss rate 
may be higher, or L36 < 0.2. 

For the wind duration , we use a M S lifetime of tMs = 8 Myr 
(Martin & Arnet3 l!995T) . iScuderi et all l fl996l) find 7V(H II) = 
4 x 10 21 cm 2 along the line of sight to Star 2, with roughly 
30% o f the extinction occuring within the LMC. IXu et alJ 
(1995) find that SN 1987 A is positioned roughly 1 kpc be- 
hind the edge of the LMC, and if one assumes Star 2 is near 
the SN, these yield an average ISM density of «o=0.4 cm" 3 . 

For L36 = 0.07-0.2, the progenitor wind will maintain an 
adiabatic shock for 0.3-1.3 Myr. As this is much shorter than 
tMs, the bubble will become radiative at a radius of 27-81 
ly. Assuming the sound speed of the undisturbed ISM is 10 
km s" 1 , the radiative shock front will reach this spee d and stall 
in ano ther ~0.8 Myr. Following the arguments in lD'Ercolel 
( 1992), adiabatic shocks become radiative at radii larger than 
those at which radiative shocks would have stalled. As such, 
the shock must recede, stalling at a new radius of 26-42 ly. 

There are subtle indications in our density measurements 
of an enhancement 80-100 ly from the SN, which may be 
the remnant of the MS-ISM snowplow. However, given the 
small expected mechanical luminosity of the progenitor, we 
question whether the Peanut is actually the remnant of that 
interaction. Alternatively, the MS bubble may be traced b y 
echoes seen roug hly 100 pc from the SN bv lXu et alJ (tl995T) . 

7.2. Red Supergiant Mass Loss 

Following the end of the MS, the gas density within its 
equilibrated bubble is expected to be quite low, and for a 
large RSG mass-loss rate, the swept-up mass within the driven 
shock will be much lower than the mass of shocked wind. The 
RSG mass-loss rate can be estimated from the mass contained 
within this structure (^1.6 T able Oil and an adopted 
RSG lifetime of 3 x 10 5 yr ( iMartin & Arnettll 19951) . yielding 
M R ~5x 10" 6 M Q yr" 1 , consistent with ISW models. A wind 
velocity of 10 km s" 1 corresponds to a mechanical luminosity 
of L36 = 1 .6 x 10~ 4 which, expanding i nto a tenuous medium, 
is classified bv lKoo & McKee dl992al) as a "slow wind" that 
becomes radiative almost immediately. Note that this rate is 
an upper limit if the MS bubble is small, as suggested above. 
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7.3. Inhomogeneities 

If the Peanut is the RSG-MS contact discontinuity, then its 
shape is a fossil record of the hydrodynamic interaction be- 
tween these two media. Here, we investigate whether the ob- 
served Peanut can be explained by a simple asymmetry in the 
RSG wind, or whether external factors need to be invoked. 
Consider the CD shell (Fig. II lb. which has an elongated, pro- 
late shape, closest to the SN just above the equatorial plane 
(r ~ 1 1 ly), and farthest toward the poles (r < 28 ly). A slow, 
radiative shock propagates as R(t)(xL x l A t l l 2 , thus for the RSG 
wind to have formed this pinched structure in the absence of 
any external inhomogeneities, the wind luminosity above the 
equato r must have been roughly 3% that at the pole. 

iD'Ercolel i 19921) find that radiative shocks stall at R m = 

22.2L^ 2 "o' /2 - For the RSG wind from Sk " 69 ° 202 t0 have 
stalled at the poles with = 1.6 x 10" 4 , the ambient MS- 
bubble density must be no = 10~ 4 cm" 3 . Similar average val- 
ues can be determined by considering the MS mass-loss rate 
and bubble size. With 3% of the polar wind luminosity, the 
equatorial wind would have stalled at 5 ly, which is less than 
half the observed position of the CD shell. It appears that the 
observed geometry of the Peanut cannot be explained by only 
an asymmetry in the RSG wind. 

One solution is that the CSM in the MS-bubble is not ho- 
mogenous, but that the MS winds had an asymmetry simi- 
lar to that of the RSG. Typically, one argues that the aver- 
age temperature inside the MS bubble should have been suf- 
ficiently high for any inhomogeneities to smooth out, mak- 
ing this scenario unlikely. An asymmetric freely-expanding 
flow into a homogenous post-shock medium will create 

an asymmtric inner termination shoc k. The simulations 

of iGarcia-Segura. Langer. & Mac Lowl i ll 9961) show that this 
shock collapses into RSG wind before the system settles anew 
into a steady flow. Thus it may be possible for a MS-wind 
asymmetry to be imprinted on subsequent outflows. An al- 
ternative is that the CD had a smaller size and asymmetry at 
the end of the RSG, and a second wind blew it into the shape 
currently observed, which we consider in greater detail in the 
next subsection. 

7.4. A Two-Loop Scenario 

Martin & Arnett ( 1995) have modelled the evolution of 
metal poor star s between 19 -21 M Q using the stellar evolu- 
tion code from lArnet3il991l) . All of these make one full blue 
loop and then return to the red, with only the models under 
20.5 Mq finishing a second loop to end as a BSG Such stars 
spend roughly 10 7 yr in the MS, 2.3 x 10 5 and 3.2 x 10 s yr 
in the first RSG and BSG, and 4.7 x 10 4 and 2.4 x 10 4 yr in 
the second RSG and BSG, respectively. We propose that the 
structures revealed by light echoes may be explained by such 
an evolutionary scenario. In particular, a long-duration BSG 
wind expanding into an asymmetric RSG environment could 
be the second wind hypothesized above. 

Qualitatively, the hourglass-shaped bubbles resulting from 
ISW simulations are poor matches to the CS hourglass, yet 
they have many similarities to the Peanut (Fig. I15> . The 
flattened, oblate (or mushroom-shaped) lobes with narrowly- 
pinched waists from simulations are very similar to the inner 
walls and conical-shell contributions to the Peanut. This is 
shown schematically in Figure^] in which the inward-facing 
shock (thin grey curve), dense contact discontinuity (thick 
gr ey curve) and vortex (g rey arrowed curve) from Figure 2 
of Martin & Arnett ( 1995) are shown at left, and compared 




FIG. 15. — Schematic showing a qualitative formation scenario for the 
Peanut. At left, the grey cur ves tra ce the favored model of BSG/RSG inter- 
action from Martin & Arnett 1 1995). Pressure or momentum directed toward 
the poles (arrows) blow this mushroom into a peanut shape (blac k dotted 
lines). At right, one quadrant of the simplified Peanut from Fig. 1131 



on the right to one quadr ant of the CD s tructure from Figure 
IT3l Of note in Martin & Arnettl (119951) . a dense vortex ap- 
pears at the outer edge of the prolate lobe, which is tempting 
to identify with the "spurs" located just outside the CD. 

We envision the following evolutionary sequence. A BSG 
accelerates wind into previously-expelled RSG material, cre- 
ating a mushroom-shaped discontinuity. This BSG wind is 
steady and long-lived, an d the interface expands self-similarly 
JPwarkadas et alJll996T) to the size of the CD/NH Peanut. 
Many mechanisms exist to direct momentum and pressure to- 
ward the polar region (indicated by black arrows in Fig. II 5I >. 
elongating the mushroom into the prolate peanut shape we ob- 
serve. These include: (1) shock focusing of the BSG wind by 
equatorial material (see Colli ns et aljU9 99). and preferential 
expansion into the lower density polar regio ns of (2) an adia- 
batic shock or (3) an ionization front (Chevalier & Dw arkadasl 
119951) . The star then loops back to the red, and for a fewx 10 4 
years fills in some fraction of the cavity blown out by the 
first BSG wind. These three mass-loss phases produce an 
extensive equatorial overdensity, which directs a substantial 
amount of the final BSG wind in the polar direction, carving 
out an hourglass-shaped nebula with no end caps. 

One major discrepancy between the ISW models and the 
CSE as revealed through echoes is that the models assume 
there is a very extensive equatorial disk of RSG material. 
Whether the equatorial material we do see is sufficient to form 
the rest of the bipolar structures must be tested in hydrody- 
namic simulations. 

In addition to providing a natural mechanism for forma- 
tion of the equatorial overdensity via a WCD, stellar rotation 
may also expl ain the progenitor's evolutionary blue loop(s) 
(Mevnet 2004), and should be revisited with hydrodynamic 
models. From Figure^] we find a maximum equator-to-pole 
density contrast of 10 between the NH equator and the fur- 
th est CD material. Th is is small compared to values found 
bv lCollins etall in"999). suggesting the evolved star may only 
have produced a wind-compressed zone and therefore did not 
have an extremely high rotation rate. However, since a rotat- 
ing star spins down as it loses mass on the MS and as it swells 
into a RSG, any significant RSG rotation requires some mech- 
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anism for spinning up the star, such as a close binary compan- 
ion. 

8. CONCLUSION 

Once it was realized that Sk -69° 202 was surrounded by 
an extensive CS nebula, it became clear that the closest ob- 
served SN in 400 years would also provide the first opportinity 
to observe the destruction of that nebula (e.g. | Luo & M cCray 
U991allBorkowski. Blondin & McCravlll997bl) . and the birth 
of a supernova remnant (SNR), both spatially resolved and in 
real-time. This event thus serves as a vital test to models of 
SN evolution, radiative shock processes, and SNR formation. 

As noted in the introduction, the ejecta-CSE interaction will 
directly illuminate the structures it impacts, and the high- 
energy emission from the resulting hot gas will illuminate 
many other regions of the CSE via photoionization. The re- 
sults from this work will be tested by this interaction, and be 
critical in interpreting the panchromatic spectacle that will en- 
sue. 

Already, the "hot spots" detailing the impa ct of t he ejecta 
with the inner edge of the ER (Lawrence et a l 12000ft confirm 
the observed spa tial offset between the SN and inner nebulae. 
Suger man et alJ J2002I) show that roughly I of the hot spots 
are located along the eastern half of the ER, while spots to- 
ward that half of the ring are also closer to the SN than those 
to the west. Both of these are suggestive of an ER that is off set 
to the west of the SN, consistent with our findings in 34.31 

An alternative interpretation to the hot spot positions and 
order of appearance is that the circumstellar medium into 
whic h the ejecta propaga te is asymmetrically denser to the 
east (iMichael et alj|20 02). Is such an imprint visible in light- 
echo data? 

Stellar winds propagating into uniform media travel faster 
and farther with increasing mechanical luminosity pv 2 
iWeaver et al.ll977t) . A higher wind density to the east would 
imply that the CD boundary separating the RSG wind from 
the MS bub ble s hould be more distant from the SN in that di- 
rection. In 34.11 we find that the CD does appear marginally 
offset from the SN toward the east. Although by no means 
conclusive, this is suggestive that an eastern asymmetry was 
present in the stellar winds as early as the RSG. 

If the RSG outflows had been denser toward the east, then 
the inner cavity which the BSG wind carved out would be 
closer to the central star in that direction as well. This would 
have resulted in an observed ER and CS hourglass that are 
offset west of the SN, and an overall density enhancement 
within the equatorial plane to the east. This qualitative argu- 
ment is again consistent with the observed geometry of the 



inner CS structures ( 34.3> . Unfortunately, the inferred densi- 
ties from 5j5]are too uncertain to permit a meaningful study of 
azimuthal asymmetries. 

In keeping with the binary-star hypotheses to explain the 
observed CSE, one could also invoke binarity as the cause 
for the east-west wind asymmetry. For example, a secondary 
in a resonant, eccentric orbit could tidally distort the pri- 
mary when near pericenter, and from the von Zeipel theorem 
(1924), the radiative flux would temporarily drop at this dis- 
tortion, with an associated local drop in wind momentum (for 
a radiatively-driven wind), yielding a longitudinal wind asym- 
metry. 

Another mechanism that could cause a longintudinally 
asymmetric wind is proper motion of the mass- losing star 
throug h the ISM. This was first trea ted by IWeaver et alJ 
i 19771) and more recen tly by e.g. Wilkin ( 199(3) and 
Comeron & Kaper ( 1998). A slow, westward motion of the 
progenitor through the ISM would pile up more material at 
the western wind-ISM contact discontinuity than to the east. 
Similarly, the star would be closer to this western discontinu- 
ity than to the west. Suc h an effect has been seen for PNe 
(Tweedv& Kwitter 1996 ) and for the LMC luminous blue 
variable SI 19 (iDanforth & Chull2001l) . A proper motion of 
less than 0.1 kms -1 is required to produce the observed of- 
fets, however, this scenario is predicated on the progenitor not 
blowing an ISM bubble during its MS that is larger than any 
of the CS structures revealed in this paper. We look forward 
to a new generation of hydrodynamical models that will read- 
dress the mechanisms of mass loss in Sk -69° 202, and sub- 
sequently the phenomena of bipolar mass loss in many more 
systems. 
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